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ABSTRACT 

The remnant of the supernova of 1006 AD, the remnant first showing evidence for the presence of X-ray 
synchrotron emission from shock-accelerated electrons, was observed for ^1000 ksec with INTEGRAL for 
the study of electron acceleration to very high energies. The aim of the observation was to characterize the 
synchrotron emission, and attempt to detect non-thermal bremsstrahlung, using the combination of IBIS and 
JEM-X spatial and spectral coverage. The source was detected with JEM-X between 2.4 and 8.4 keV bands, 
and not detected with either ISGRI or SPI above 20 keV. The ISGRI upper limit is about a factor of four above 
current model predictions, but confirms the presence of steepening in the power-law extrapolated from lower 
energies (< 4 keV). 

Subject headings: ISM:individual (SN1006), supernova remnants, X-rays:observations, radiation 
mechanismsinon-thermal 



1. INTRODUCTION 

Supernova remnants (SNRs) have long been thought to be 
the primary site of Galactic cosmic ray acceleration up to 
the "knee" feature in the integrated cosmic -ray spectrum near 
3000 TeV, as the supernova shocks are one of the few mech- 
anisms that could provide enough energy to support this pop- 
ulation ( iDveret al.l2001h . However, many features of the ac- 
celeration process, including injection physics, efficiency, and 
maximum electron and ion energies, are not yet clear Hard 
X-ray observations (above 10 keV) of SNRs may cast light on 
this poorly understood process. 

High energy electrons in SNRs produce X-rays via two 
mechanisms, non-thermal 

bremsstrahlung and synchrotron radiation. Electrons can also 
produce gamma-rays up to TeV energies by inverse Comp- 
ton (IC) scattering of any photons present, such as the cos- 
mic microwave-background (CMB). In addition, relativis- 
tic protons can produce gamma-rays from the decay of tt" 
particles from inelastic ion-ion collision. All these pro- 
cesses have been extensively modeled for SNRs by vari- 
ous g roups (Sturner et al. 1997; Gaisser, Protheroe & Stane^ 
IBaring 



ll998HBaring et al...l999.) . In hard X-rays, synchrotron radi- 
ation from the tail of the electron distribution may compete 
with non-thermal bremsstrahlung from the very lowest-energy 
accelerated electrons. 

SN 1006 has been the prototype laboratory for the study of 
electron acceleration to high energies in shocks. X-rays from 
this object were first reported by Winkler & Laird ( 1976). The 
earlier featureless spectrum (Becker etal. 1980) was mod- 
eled as the loss-steepened extrapolation of t he ra dio syn- 
chrotron spectrum by iRevnolds & Chevalied (Il98ll) . Later, 

' Space Sciences Laboratory, 7 Gauss Way, University of California, 
Berkeley, CA, 94720-7450, USA. 

^ Sabanci University, Orhanli-Tuzla 34956, Istanbul, Turkey. 
Department of Physics, NC State University, 2700 Stinson Drive, Box 
8202, Raleigh, NC 27695, USA. 

* Department of Physics, University of California, 366 Le Conte Hall, 
Berkeley, CA, 94720-7300, USA. 

' Danish National Space Center, Juliane Maries Vej 30, DK-2100 Copen- 
hagen 0, Denmark. 

^ Service d'Astrophysique, CEA-Saclay, 91191, Gif-Sur-Yvette, France 

' APC-UMR 7164, 11 place M. Berthelot, 75231 Paris, France 

^ SIRTF Science Center, California Institute of Technology Mail Stop 220- 
6, Pasadena, CA 91125 



observations by ASCA ("Kovama et al."'1995") showed that the 
limbs have featureless spectra well described by power-laws, 
whereas the interior has a thermal, line-dominate d spectrum. 
The source was also observed with RXTE, and iDver et all 
(2001) showed that e laborate synchrotron emission models 
(Reynolds 1996, 1998) fit the combined RXTE - ASCA spec- 
trum reasonably well. 

Electrons producing keV synchrotron emission could also 
produce very high-energy photons (in the TeV range) by IC 
upscattering of CMB photons (Pohl 1996). The TeV flux de- 
pends on the electron distribution, and in conjunction with the 
synchrotron flux, a mean magnetic field strength of the rem- 
nant can be deduced. A detection of the northeast (NE) limb 
of SN 1006 was reported in ground-based TeV observations 
by CAA^GA/?(9(9-/ (Tanimori et al. 1998). It is interesting that 
only one limb was detected, even though the X-ray spectra 
of the two limbs are similar (Allen. Petre & Gotthelf 200 1^ 
Dver, Revnolds & Borkowski 2004). To have a discrepancy 
in the TeV band, the electron spectra, magnetic field strengths, 
or synchrotron and IC emission-an gle distributions of the 
two rims would have to be different (* A llen. Petre & Gotthelfl 
1^01), but in a way that does not produce significant differ- 
ences in the X-ray band. 

The nature of X-ray emission from SN 1006 above 10 keV 
is still uncertain. Below 10 keV, synchrotron emission is the 
most plausible explanation. For synchrotron radiation, the 
quantitative inferences apply only to the exponential cutoff of 
the electron distribution. On the other hand, most of the ac- 
celerated electrons and much of their total energy reside in the 
lowest energy non-thermal electrons, whose bremsstrahlung 
emission could become dominant above 30 keV. In princi- 
ple, the INTEGRAL observatory can examine the effects of 
both the lowest and highest energy non-thermal electrons by 
distinguishing the synchrotron and bremsstrahlung emission 
w ith its imaging an d spectral capabilities. 

lRevnoldsl ( ll999 t) modeled bremsstrahlung and synchrotron 
emission from SN 1006. Synchrotron hard X-rays should 
be concentrated in two bright opposing limbs like the ra- 
dio emission, and should dominate the emis sion below 30 
keV The images and spectra taken wi th ASCA dKovama et alJ 
Chandra X-ray Observatory ("Lon g et aLl l200?). and 
XMM -Newton (|Rpthenflug et al. 2004) confirm this below 10 
keV. The non-thermal bremsstrahlung, resulting from slightly 
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supra-thermal shock-accelerated electrons interacting with 
thermal ions, is likely to be more symmetrically distributed, 
and will dominate at some energy between 30 and 300 keV. 
The bremsstrahlung flux will scale with the product of the 
thermal gas density n± and the relativistic-electron density 
Wciei- The former can be constrained by observations of ther- 
mal X-ray emission, while the latter can be deduced from ra- 
dio synchrotron fluxes if the magnetic field is known. 

The TeV spectrum reported from CANGAROO-I 
(Ijanimori et al. 1998) could be well-described by IC 
upscattered CMB photons, using a power-law electro n spec- 
trum w ith an exponential cutoff as described in Dv er et al.l 
This fit gave a post-shock magnetic field of about 
10 ^G. However, this result is now called into question by 
the observations of H.E.S.S., which did not detect the source 
despite better sensitivity compared to CANGAROO-I. The 
TeV upper limits from H.E.S.S. ar e about a factor of 1 
below the CANGAROO-I results'* (lAharonian et all 120051) . 
These limits constrain IC upscattering of CMB photons by 
the same electrons that produce X-ray synchrotron emission: 
tighter limits mean fewer electrons, a higher magnetic field 
(Aharonian et al. place a lower limit of 25 fiG on the post- 
shock magnetic field), and less non-thermal bremsstrahlung. 
Lowering the possible IC-CMB flux by a factor of 10 directly 
lowers the allowable relativistic-electron density by the same 
factor, and hence lowers the predicted bremsstrahlung flux by 
an order of magnitude. 

Even if bremsstrahlung is not detected, the detailed shape of 
the steepening synchrotron spectrum can provide information 
on the physical process causing the cutoff in the electron spec- 
trum. This is crucial information for understanding the accel- 
eration of cosmic rays, since if the cutoff is due to radiative 
losses on electrons, the proton spectrum might extend to much 
higher energies, perhaps as high as the knee. However, if the 
finite remnant age (or size) or some change in diffusive prop- 
erties of the upstream medium causes the cutoff, it should also 
cut the proto n spectrum off at a similar energy (between 10 
and 100 TeV: lDver et al.ll2001 ). far below the "knee" energy. 
Detailed models show subtle but potentially distinguishable 
differences in the shape of the synchrotron spectrum above 10 
keV, with loss-limited spectra (due to higher magnetic fields) 
being somewhat harder. 

Our group has observed SN 100 6 for ^10 00 ks with the 
INTEGRAL Observatory ( Winkler eTal]l2003l) in AO-1 with 
the main aim of detecting and characterizing synchrotron 
emission, and distinguishing synchrotron and non-thermal 
bremsstrahlung emission by comparing the IBIS/ISGRI and 
JEM-X images to the model images. In this work, we will 
discuss the results of the analysis of the INTEGRAL data, and 
place limits on the synchrotron and bremsstrahlung emission 
from SN 1006. 

2. OBSERVATIONS AND ANALYSIS 

The INTEGRAL observations took place in two sets. The 
^250 ks first set ("Set I") was conducted early in the mission, 
between Jan, 11, 2003, and Jan, 20, 2003, corresponding to 
INTEGRAL revolutions 30 and 32. The '^750 ks second set 
("Set 11") was conducted between Jan, 20, 2004 and Jan, 30, 
2004, during revolutions 155-158. These two sets have differ- 
ent observational characteristics for different instruments as 
explained below, and only the data from Set II are used for 

' Recent measurements by CA NGAROO-UI has also claimed null result 
on SN 1006 ITanimori et all2005l) . 



this paper. We did not use SPI (Spectr omet er on INTEGRAL, 
Vedrenne et al. 2003), as ISGRI (see § l2.2l for more informa- 
tion) places much stricter limits in the hard X-ray band. Be- 
fore the general analysis for all instruments, we filtered out 
the pointings with high Anti-Coincidence Shield rates, mostly 
occurring during the entry and exit of the radiation belts. 

2.1. The JEM-X analysis 

The Joint European X-ray Monitor, JEM-X, consists of two 
identical high pressure imaging micro-strip gas chambers, and 
makes observations simultaneously with the main instruments 
on INTEGRAL, albeit with a narrower fully-coded field of 
view of 4.8°. The energy band is 3- 35 keV and the angular 
resolution is 3.35 ' JLund et alJl2003l) . Due to a problem with 
eroding anodes, the high voltage in the JEM-X detectors was 
lowered, and a new background rejection criterion was imple- 
mented after the launch. One of the detector pair is being kept 
in a safe state, and during our observations only JEM-X 2 was 
operational. 

Set I was conducted before the new background rejection 
criteria were implemented in JEM-X, and therefore was not 
included in this analysis. We note that, although the total ex- 
posure for Set II is ~750 ks, the effective exposure time of the 
central object is approximately 250 ks due to the vignetting 
of the JEM-X instrument during the 25 point dither We have 
obtained JEM-X images in 4 energy bands using the JEM- 
X Midisky offline software package available from the DNSC 
(Lundetal. 2004). These energy bands"' are 2.4-4.2 keV, 
4.2-8.4 keV, 8.4-14 keV, and 14-35 keV. The images from 
each pointi ng are then mosaic ked using the mosaic-weight 
program ( Chenevez et al. 2004). 

The fluxes and the significance values shown in Tableware 
derived using the inner ASCA contours enclosing the NE and 
SW limbs to define two shape templates. The NE template 
is constructed with 23 pixels inside the inner ASCA contour 
seen in Fig.^on the left. The SW template is also constructed 
similarly, using the 19 pixels inside the inner ASCA contour 
on the right. The count excesses in the JEM-X mosaic im- 
ages are then determined inside these two regions. The rel- 
evant noise figures are derived by defining a number of non- 
overlapping regions with the reference templates within a 60 
X 60 pixel field (90 x 90 arcminutes) centered on SN1006. 
114 NE-templates, each containing 23 pixels, can be arranged 
within the 60 x 60 field. For the S W-template the correspond- 
ing numbers are 1 12 regions, each with 19 pixels. The statisti- 
cal properties of the excess counts in these two sets of regions 
are used to derive RMS-noise of the background for regions 
of the two shapes. The signal-to-noise ratio derived in this 
way will obey normal statistics because the region templates 
are defined independently of the JEM-X data. 

The fluxes are derived by comparing the excesses in the 
SN 1006 mosaics with corresponding excesses in mosaic im- 
ages of the Crab Nebula obtained with the same INTEGRAL 
dither pattern, and correcting for the difference in the effective 
observation time. 

2.2. The ISGRI analysis 

One of the two main instruments on INT EGRAL, IBIS (Im- 
ager on Board the INTEGRAL Satellite, Ubertini et al. 200?), 
consists of two cameras. The INTEGRAL Soft Gamma-ray 

The low energy threshold is yariable over different parts of the JEM-X 
detector. 
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Imager, ISGRI, is the low-energy camera of the IBIS tele- 
scope (Lebrunetal. 2003). It has a large sensitive area of 
2621 cm^ made up of 16384 CdTe pixels. The angular reso- 
lution is ^13' (Gros et al. 2003), and the fully-coded field of 
view is 9°. The energy range is 20 keV - 10 MeV. The older 
background maps for Set I observations resulted in much 
noisier images compared to Set II images, and therefore we 
limit the analysis to the data from Set II. We used OSA 4.2 
(jpoldwurmet al. 2003) standard programs to obtain images 
in 20-40 keV and 40-100 keV bands. 

Since the ISGRI system point spread function is ~13', 
SN 1006 ('^30' diameter) appears as an extended object to 
the imager Estimation of the flux of such sources with a 
coded-mask instrument is complicated, since the mask pat- 
terns used in gamma-ray astronomy are optimized for point 
sources. One needs to use simulations to obtain the effect of 
the extended nature of the source on the image and the mea- 
sured flux. Such simulations have been conducted for ISGRI 
for different extended source geom etries, including SN 1006 
(iRenaud et aVlOOS). For SN 1006, Renaud et al. (2005) used 
an input image of the expected synchrotron emission map 
in 20-40 keV band based on the simulations described in 
iRevnoldsl fl999,) . The principle of the ISGRI simulations is 
the following: For each point-like source constituting the ex- 
tended one, the corresponding shadowgram (the image of the 
mask pattern illuminated and projected onto ISGRI) is calcu- 
lated. The final expected shadowgram of SN 1006 is obtained 
by summing all these contributions. Then, the standard de- 
convolution (see Goldwurmet al. 2003, for details) in OSA 
is applied to obtain the reconstructed image. With this tech- 
nique, Renaud et al. (2005) obtained a reduction factor of 0.7: 
that is, ISGRI would detect 70% of the true flux at each limb. 

3. RESULTS 

3.1. JEM-X results 

The source is detected at the limbs in 2.4-4.2 keV and 
4.2-8.4 keV bands (see Fig. [T] and Table Q. This is the 
first time that the structure of an individual extended 
source has been imaged with INTEGRAL. It appears that 
the South- West (SW) limb (the limb on the right in the im- 
ages) is stronger in the 2.4—4.2 keV band at about the la 
level, but the trend is reversed at higher energies. At 4.2- 
8.4 keV, the NE limb is stronger, a result with higher signif- 
icance. An excess at the position of the NE limb is present 
in the 8.4-14 keV band (1.7 cr), but no excess is seen in the 
SW limb. We note that XMM-Newton data also points to an 
asymmetry in flux coming from the NE and SW limbs such 
that the NE limb^ets r elatively stronger as energy increases 
(|Rothenflu2 et al. 2004^ A similar trend is seen in ASCA data 
(pDver. Revnolds & Borkowski 2004), but the differences are 
small in the ASCA band (below 8 keV). The small excess of 
the SW over NE limbs in the 2.4^.2 keV band we see is not 
supported by ASCA or XMM-Newton observations in that en- 
ergy range. 

Figure |2] shows our combined JEM-X fluxes for both 
limbs, compared with previous sp atially in te grated 
pre-A^CA flux measurements (see iRevnoldi Il996t 
iHamiltonTs arazin & Szvmko wiakI 119861 for references). 
The integrated fluxes have been converted to flux densities 
assuming a photon index F of 3.0 (Allen. Petre & Gotthelf 
|2P01). Also shown are a pair of model curves for the escape 
model that provided a good fit to A SCA and RXTE data in 
iDver. Revnolds & BorkowskH (120041) . The parameters of the 




Fig. 1 . — JEM-X reconstructed images of SN 1006 in three energy bands. 
Top: 2.4^.2 keV, middle: 4.2-8.4 keV band, bottom: 8.4-14 keV. The ASCA 
contours are overlaid. The detection significances, fluxes, and upper limit 
fluxes are shown in TablelTI (The electronic version has this figure in color.) 



best fit are a roll-off frequency of Vioii = 3.0 x 10 Hz and 
an electron energy index of 2.2 (implying a radio spectral 
index a = 0.6 or radio photon index F = a + 1 = 1.6). These 
parameters imply an e-folding energy of the exponential 
cutoff (due to escape) of 32(B2/10/iGauss)"'/^ TeV (where 
Z?2 is the post-shock field). The curves in Figure|2]correspond 
to ±1(7 errors on i/roii: (2.8-3.1) x 10'^ Hz. The dashed 
line indicates the bremsstrahlung prediction for an upstream 
density of 0.2 cm"-' and B2 = lO/iGauss. 

3.2. ISGRI results 

We obtained ISGRI images in different energy bands using 
OSA 4.2. SN 1006 was not detected in any band. The 20- 
40 keV sigma image is shown in Fig.|5] The 3 cr upper limit 
(sensitivity limit) for a point source for ^750 ks observing 
time is ^9 x 10"^ photon cm"^ s"'. If synchrotron dominates 



TABLE 1 
JEM-X Summary 



North-East Limb 








Energy band 


Flux 


Statistical 


Model flux" 


(keV) 


(10^* photon/cm- s) 


significance 


(lO^** photon/cm- s) 


2.4^.2'' 


10 ±4 


2.6 


14 


4.2-8.4 


11 ±2 


5.0 


5.9 


8.4-14 


3.3 ± 1.9 


1.7 


1.4 


South- West Limb 








2.4^.2 


15 ±4 


4.0 


13 


4.2-8.4 


3.6 ±2.3 


1.6 


4.8 


8.4-14 


<5.7 


3^ 





"Flux from extrapolated model fi t to the ASCA spectrum (0.8-9 keV), 
IDver. Reynolds & B orkowski'2004 

*The actual low energy threshold is variable over different parts of the JEM- 
X detector. 

' 3(T upper limit 
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Fig. 2. — Integrated hard X-ray spectrum of SN 1006. The heavy bars show 
the energy bins over which our JEM-X fluxes are integrated. The integrated 
fluxes were converted to spectral fluxes assuming a power-law index of 3.0, 
and plotted at the median energy of each bin. The errors in the bin energies 
are insignificant (much smaller than the line widths). The upper limit bar 
is from ISGRI. Open circles are pre-A5CA (see text) observations, and the 
two solid lines are the ± Icr model fits to ASCA and RXTE data, extrapolated 
to the higher energies (Dyer et al. 2001). The dotted-and-dashed line is the 
extrapolation from the Chandra spectral fit (Loneetal. 2003). The dashed 
line is the bremsstrahlung prediction of the model described in text. 



in 20^0 keV band, the emission would be concentrated in 
two limbs. The extended natur e of the source causes a reduc- 
tion factor of 0.7 in flux (see § l2.2t . Therefore the 3 a upper 
limit for a synchrotron dominated source is 1.3 x 10""* pho- 
ton cm"^ s"' at each limb. For synchrotron dominated emis- 
sion, the expected total flux in the 20-40 keV band (based 
on the models shown in Figure^ is (2.9-3.6) x 10"^ photon 
cm"^ s~' for the entire remnant, or roughly half this at each 
limb. The bremsstrahlung emission from the model shown in 
Figure 12] using the lower- limit 25 fiG downstream magnetic 
field from H.E.S.S. TeV observations, is predicted to be about 



2 X 10 photon cm s" , far below ISGRI's sensitivity. In 
fact, in this case, since the emission would be coming from 
a larger area, the upper limit set by our observations is even 
higher. The predicted crossover energy where bremsstrahlung 
and synchrotron emission become comparable is in the vicin- 
ity of 200 keV. 

Long et al. (2003) reported that a power-law photon index 
of 2.30 fit the Chandra data well between 0.5 and 5 keV for 
the emission at the limbs. An extrapolation of this fit from 5 




Fig. 3.— The ISGRI sigma image around SN 1006 in 20-40 keV band. 
SN 1006 is not detected. A nearby source (possibly VV 780, 11.6 cr detec- 
tion lKalemci. B oggs & Lund 2005) is also shown for comparison. The A5CA 
contours are overlaid to show the expected position of SN 1006. The feature 
at the center of SN 1006 is not significant, and is possibly due to inaccurate 
background map. (The electronic version has this figure in color.) 



keV to the ISGRI range is also shown in Fig.|2]with a dotted- 
and-dashed line. Our ISGRI upper limit is about a factor of 2 
below the extrapolation of this power-law to 28 keV. 

4. DISCUSSION 

Coded-mask imaging is a complex process, and optimal 
techniques are still being developed to clean noisy images, 
both for JEM-X and ISGRI. Our results confirm and extend 
results obtained with A5CA, Chandra X-ray Observatory, and 
XMM-Newton, and show that relatively small increases in sen- 
sitivity may allow ISGRI to detect predicted synchrotron ra- 
diation for some models. 

The prospect of detecting bremsstrahlung X-rays or 
gamma-rays from SN 1006 has become considerably more 
remote with the combination of lower estimates for the ambi- 
ent density and the much higher magnetic field that would 
be required to explain the lack of IC from CMB photons 
in the H.E.S.S. observations. However, it is quite possible 
that a somewhat longer INTEGRAL observation scheduled for 
Cycle 3 will have sufficient sensitivity to detect synchrotron 
emission in the 18-40 keV band. A model that can describe 
the observations of Figure |2l but with a strong enough mag- 
netic field to satisfy the H.E.S.S. constraints (downstream 
B > 25/iG) is somewhat harder than the spectra shown in Fig- 
ure|2l and predicts a flux in the 18-40 keV band of 0.3 x 10""* 
photon cm"^ s"'. Some models predict considerably lower 
fluxes in the 18-40 keV band, so ISGRI detection would not 
only extend and confirm the presence of hard synchrotron X- 
rays, but could provide useful model discrimination. 

We have obtained the first observations with JEM-X of an 
extended source, between 2.4 and 14 keV. The fluxes we de- 
rive are consistent with those of earlier imaging observations 
where they overlap, and support the identification of the con- 
tinuum emission of SN 1006 as synchrotron radiation from a 
slowly dropping off electron distribution. Our ISGRI 3a up- 
per limit is about a factor of 4 higher than the prediction of the 
model that best fits the soft X-ray continuum. With ISGRI, we 
confirm at higher energies than has been previously reported 
that the X-ray spectrum of SN 1006 must be steepening. A 
somewhat longer observation, and more developed data anal- 
ysis techniques, should allow the detection of SN 1006 in the 
18-40 keV band with ISGRI, and above 8.4 keV with JEM-X, 
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and can provide important modeling constraints. 
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